Heisenberg antiferromagnet with anisotropic exchange on the kagome lattice: 
Description of the magnetic properties of volborthite 
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We study the properties of the Heisenberg antiferromagnet with spatially anisotropic nearest- 
neighbour exchange couplings on the kagome net, i.e., with coupling J in one lattice direction and 
couplings J' along the other two directions. For J /J' > 1, this model is believed to describe the 
magnetic properties of the mineral volborthite. In the classical limit, it exhibits two kinds of ground 
states: a ferrimagnetic state for J/ J' < 1/2 and a large manifold of canted spin states for J/ J' > 1/2. 
To include quantum effects self-consistently, we investigate the Sp(JV) symmetric generalisation of 
the original SU(2) symmetric model in the large-TV" limit. In addition to the dependence on the 
anisotropy, the Sp(7V) symmetric model depends on a parameter k that measures the importance of 
quantum effects. Our numerical calculations reveal that in the k-J/J' plane, the system shows a rich 
phase diagram containing a ferrimagnetic phase, an incommensurate phase, and a decoupled chain 
phase, the latter two with short- and long-range order. We corroborate these results by showing that 
the boundaries between the various phases and several other features of the Sp(jV) phase diagram 
can be determined by analytical calculations. Finally, the application of a block-spin perturbation 
expansion to the trimerised version of the original spin- 1/2 model leads us to suggest that in the 
limit of strong anisotropy, J/ J' ^> 1, the ground state of the original model is a collinearly ordered 
antiferromagnet, which is separated from the incommensurate state by a quantum phase transition. 

PACS numbers: 75.10. Jm,75.30.Kz, 75. 50.Ee 
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I. INTRODUCTION 

In the ongoing search for novel states of condensed mat- 
ter, frustrated antiferromagnets have played a key role 
(for a recent review, see Rcf. 1). Among the many sub- 
stances that have been investigated experimentally and 
the numerous spin models that have been studied theo- 
retically, those in which the magnetic ions occupy the ver- 
tices of corner-sharing frustrating entities have attracted 
particular attention in this context. The best known ex- 
amples are the kagome antiferromagnet (KAF), consist- 
ing of corner sharing triangles, and the pyrochlore anti- 
ferromagnet, consisting of corner sharing tetrahedra (see 
Fig. 1). 

The main distinction between the KAF, the pyrochlore 
antiferromagnet and other frustrated and unfrustratcd 




(a) (b) 
FIG. 1: kagome lattice (a), pyrochlore lattice (b) 



magnets is the large ground-state degeneracy of the for- 
mer: classical Heisenberg antiferromagnets with nearest- 
neighbour interactions on corner-sharing lattices have a 
large ground-state degeneracy, which in the above two ex- 
amples even leads to a finite ground-state entropy (see, 
e.g., Ref. 2 and references therein). Quantum effects 
may lift this degeneracy, and, indeed, in numerical stud- 
ies of small cells of the spin | KAF, an exponentially 
large number of very low-lying quantum states has been 
observed 3,4 . It has been suggested that this abundance 
of low-lying states can be understood in a description of 
the low-energy physics of the quantum KAF as spin liquid 
consisting of nearest-neighbour spin singlets 5,6 . However, 
a complete picture of the ground state and of the excita- 
tions of the KAF is still missing. Further theoretical, but 
also experimental studies with emphasis on the quantum 
properties of the KAF arc therefore highly desirable. In 
this last respect, the mineral volborthite is a very promis- 
ing candidate. It has been the subject of several recent 
experimental investigations 7-10 . The magnetic lattice of 
this natural antiferromagnet consists of the 5 = 1/2 spins 
of Cu 2+ ions that are located on the vertices of well sep- 
arated planar kagome- like nets. A monoclinic distortion 
of the lattice leads to a slight difference between the ex- 
change couplings along one lattice direction (J) and the 
two other directions (J') (sec Fig. 2). Since neither signs 
of long-range order nor signs of a spin-gapped singlet 
ground-state were found in experiments on volborthite, 
the substance seems to be a good candidate for the ob- 
servation of the low-energy features that are thought to 
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FIG. 2: (Color online) Anisotropic kagome model. The cou- 
pling J' and the nearest neighbour distance will be set equal 
to unity in the calculations. 6i(\\ e x ), 82, and £3 are the three 
primitive lattice vectors of the kagome net. 



be typical for kagome type antiferromagnets 1 . 
Whether and to what extent the different exchange cou- 
plings along different lattice directions of the kagome net 
of volborthite influence the low-energy physics of the sys- 
tem is presently unknown. In the present paper, we study 
this question on the basis of the model Hamiltonian 



(1) 



(k,i) 



The symbols and (k,i) denote, respectively, bonds 
between nearest-neighbour sites on the horizontal chains 
(a, b) and bonds between the middle sites (c) and the 
sites a, b, see Fig. 2. Since the physics of this model de- 
pends only on the ratio J/J' of the exchange constants, 
we set J' = 1 in the sequel. We will consider the spatially 
anisotropic kagome antiferromagnet (AKAF), Eq. (1), in 
the full range of J, < J < 00 since this is of theoretical 
interest: one expects to see quantum phase transitions 
as J is increased. It is of particular interest to find out 
whether there is a transition from two-dimensional mag- 
netic states to a set of decoupled chains with free spins 
on the axes between the chains for large values of J. 
The paper is organised as follows. In Sec. II, we consider 
the model (1) in the classical limit. At this level, we find 
no sign of a transition from the two-dimensional magnet 
to a set of decoupled chains as J increases to infinity. 
Nonetheless, the ground-state degeneracy, as well as the 
spin wave spectrum are found to change qualitatively 
as the anisotropy of the model varies. In Sec. Ill, we 
consider a generalisation of the SU(2) symmetric model 
(1) to the Sp(jV) symmetric version 11 ' 12 and describe 
its properties in the large- M limit, where a mean-field 
treatment of the model is adequate. We obtain a de- 
tailed description of how possible ground states of the 
model depend on the coupling J and on the spin length 
S. A fairly rich phase diagram with a ferrimagnetic 
phase for small J, long-range ordered and short-ranged 
incommensurate phases for intermediate values of J, and 
a decouplcd-chain phase for large J emerges. Parts of 
these results have been published previously, sec Ref. 13. 



In Sec. IV, we devise trial quantum ground states of the 
original S — 1/2 model. We chose the states such that 
they are exact eigenstates of Hakaf , if the couplings on 
the upward pointing triangles of Fig. 2 are switched off, 
and we then treat these couplings perturbatively. In the 
limit J — > 00 this yields an effective Hamiltonian for 
the spins on the c sites which represents an anisotropic 
triangular antiferromagnet. The conclusions of Starykh 
and Barents 14 about the ground state of this effective 
model lead us to conjecture the existence of a quantum 
phase transition in the AKAF for large J. In Sec. V, we 
summarise and discuss our results. In two Appendices, 
we present technical details of the counting procedure for 
the classical ground-states, and of the Ginzburg-Landau 
type procedure that allows us to determine the bound- 
aries in the phase diagram analytically. 



II. CLASSICAL AND SEMICLASSICAL 
ASPECTS 

Similar to other isotropic spin models on lattices with 
triangular elementary cells, the classical ground states of 
Wakafj Eq. (1), are spin configurations, which satisfy the 
condition that for each elementary triangular plaquette 
of the lattice, Fig. 2, the energy is minimal. 
For J = 0, this yields a ferrimagnetic state with the chain 
spins aligned in one direction and the middle spins point- 
ing in the opposite direction, so that the total magneti- 
sation is M = N^yS (-/V v : number of downward pointing 
triangles, iV v = 7V s /3 where N s is the number of sites of 
the system). We illustrate this situation in Fig. 3. Ac- 
cording to the Licb-Mattis theorem, the exact quantum 
ground state (GS) of the model 7Yakaf also has total 
spin S tot = N^S for J = 0, see Ref. 15, i.e., for J = 0, 
the quantum GS is ferrimagnetic too. By continuity, one 
expects the quantum GS to remain ferrimagnetic for suf- 
ficiently small finite J. This will be confirmed by our 
considerations of the large- limit of the Sp(AT) version 
of our model (see the analytical and numerical work in 
Sects. Ill, III C and Appendix B) and by the block spin 
perturbation approach (Sec. IV). Classically, the ferri- 
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FIG. 3: (Color online) Ferrimagnetic state for J = 0, 
i.e., when there is no coupling between chain spins, cf. Fig. 2. 
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FIG. 4: (Color online) Spin- wave frequencies w a (q), a — 3 

1,2,3 for J = 0.4; the contour at the top of the plot marks FI q. 6: (Color online) Same as Fig.4 for J = 0.6. 

half the Brillouin zone. 



magnetic state remains stable up to J = 1/2. The ex- 
citation spectrum of the ferrimagnetic state obtained in 
linear spin- wave (LSW) approximation is shown in Fig. 4. 
The analytic expressions for these three frequency sur- 
faces w a (q), a = 1,2,3, are obtained as solutions of a 
third-order secular equation and are too lengthy to be 
presented here. However, one can easily assure oneself 
that the dispersion of the gapless mode is quadratic at 
the origin. Thus, one has the typical mode structure of 
a fcrrimagnet here with one ferromagnetic mode and two 
optical modes, see, e.g., Ref. 16. As J increases towards 
1/2, the ferromagnetic frequency surface looses its disper- 
sion and turns into a plane of zero modes, one zero mode 
for each wave vector in the magnetic Brillouin zone (BZ), 
at J = 1/2. The gap of the lower optical mode closes at 
this value of J in the centre of the BZ and the dispersion 
of this mode becomes linear for small wave vectors as for 
an antiferromagnetic spin-wave mode. 
At J = 1/2, the classical GS configuration changes from 
the unique ferrimagnetic state to an ensemble of degener- 
ate canted coplanar states. These states are characterised 
by two variables: the angle 8, which the middle spin of 
a given triangular plaquctte forms with the two chain 
spins of the same plaquette (see Fig. 5), and the two 
valued chirality x = ilj which denotes the direction in 
which the spins turn as one moves around the plaquette 
in the mathematically positive sense. 
For J > 1/2, the requirement that the energy of any of 
the elementary triangular plaquettes of the lattice Fig. 2 
be minimal is 6 = arccos(— 1/(2 J)), (9 > 0). The differ- 
ent degenerate canted states arise from different possibil- 
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FIG. 5: (Color online) Canted spins of the AKAF at J > 1/2. 



ities to assign positive or negative chiralities to the pla- 
quettes of the lattice. We show in the Appendix A that 
for the general case of 8 ^ 27r/3 ( J ^ 1), the number of 
spin configurations, iVg™ so does not grow exponentially 

with the number of sites. Rather, N£ n s ls ° < 2 a V^, 
where a < 3. This implies that the ground-state en- 
tropy per spin of the classical AKAF vanishes in the 
thermodynamic limit. In this respect, the anisotropic 
model differs qualitatively from the isotropic KAF in 
the classical limit, which has an extensive entropy per 
spin. In the limit J — > 1, the anisotropic model ap- 
proaches the isotropic KAF. Hence, one expects that for 
the anisotropic model there is an extensive number of 
low-lying excited states that become degenerate with the 
GS in the isotropic limit. 

As in the case of the isotropic KAF, the spin- wave Hamil- 
tonian is in linear order independent of the particular 
classical GS that has been chosen as the starting point 
of the expansion, Ref. 17. This implies that lowest-order 
quantum fluctuation do not select one or a group of clas- 
sical GSs as true GSs, i.e., the possible ordering effects 
of quantum fluctuations are not captured by the linear 
spin-wave (LSW) approximation. Figs. 6, 7 show the 
spin- wave frequency surfaces for J = 0.6 and for J = 3. It 
is easy to show analytically that, as is illustrated in these 
figures, the plane of zero frequency modes persists for all 
values of J greater than 1/2. The surfaces for J < 1/2 
and for J > 1/2 join smoothly at J = 1/2. Thus, in the 




FIG. 7: (Color online) Same as Fig.4 for J = 3. 



4 



LSW approximation, the transition from the ferrimag- 
netically ordered state to the canted spin states appears 
to be of second order. For J ^> 1, the nonzero frequen- 
cies gradually loose their dispersion perpendicular to the 
strong- J direction and take the shape of the spin-wave 
spectrum of antiferromagnetic chains parallel to this di- 
rection. However, no sign of a further transition from 
the canted spin states to a set of decoupled spin chains is 
found in this semiclassical approach. In the next section, 
we will consider the symplectic Sp(A/") generalisation of 
the antiferromagnetic model 7Yakaf in the large- A/" limit. 
This approach, which was first proposed by Read and 
Sachdev, Refs. 11,12, as a method to study frustrated 
antifcrromagnets, has the benefit of including the or- 
dering effects of quantum fluctuations self-consistently. 
It is of particular interest for spin models with two or 
more competing exchange couplings in the different lat- 
tice directions or over different lattice distances such as 
the present model, the J1-J2-J3 model 11 , the Shastry- 
Sutherland antiferromagnet 18 and the anisotropic trian- 
gular antiferromagnet 19 . For these models, it has pro- 
vided an unbiased selection of possible GSs that may or 
may not be ordered depending on the value of a param- 
eter n, which is connected with the spin length S (see 
below). 

III. MEAN FIELD Sp(A0 APPROACH 

A. Brief review of the method 

For a general antiferromagnetic Heisenberg model with a 
positive interaction matrix Jy, 

H = Jij Si ■ Sj , (2) 

i>j 

the Sp(AT) generalisation reads 

n sm) = - £ • (3) 

i>j 

Here. 



J = £ (4) 



is the 2Af x 2Af generalisation of the 2x2 antisymmetric 
tensor 

~(-°i + o). < 5 » 



and b" with a = 1, . . . , 2Af are the Sp(7V) boson anni- 
hilation operators. (Here and in the sequel, we closely 
follow the notation of Ref. 12; in particular, summation 
over repeated upper and lower indices is implied.) Thus, 
J a ^b\ a bjp is the generalisation of the creation operator 

£ al3 bl a bjp for a singlet on the bond (i, j). For the special 
case Af = 1 , one finds 

(J a0 bt a b] )(J 1 sb]b 5 j ) = -2Si ■ Sj + n H n bj /2 + 5 t] n u , 

(6) 

where 

nu = bib? (7) 
is the boson number operator at site i and where 

S< = (8) 

is the usual SU(2) spin operator at site i. (r are the 
Pauli matrices). Then, if one imposes the constraint 
that the number of bosons is the same for all lattice sites, 
= rib, the Hamiltonian Wgp(i) is the familiar SU(2) in- 
variant antiferromagnetic Heisenberg Hamiltonian (plus 
some constants) with rib = 2S. 

In the subsequent exposition, we shall consider a Hamil- 
tonian of the form (3) in the large- AT limit following the 
strategy of Refs. 11,12. Depending on the values of the 
couplings and of k, the GS of 7isp(A0 ma Y either 
break the global Sp(W) symmetry and exhibit LRO or 
it may be Sp(jV) symmetric with only SRO. Breaking of 
the Sp(jV) symmetry will happen through condensation, 
i.e., by macroscopic occupation of one of the Bose fields 
b a . To allow for this, we introduce the parametrisation 




with a = (ma), m = 1, • • • , Af, m = 2, • • • , AT and 
a =t, I- The field x? is proportional to the condensate 
amplitude, = V^AA 6™ xf. Aiming at a mean field 

treatment of the Hamiltonian 7isp(AA)i which becomes ex- 
act in the large Af limit, we decouple the quartic part by 
the Hubbard-Stratonovich technique with complex fields 
Qij = —Qji an d with Lagrange multipliers that en- 
force the local constraints (7). The variables Qij which 
are defined on nearest neighbour bonds of the lattice are 
expectation values of the bond singlet creation opera- 
tors in the GS, Q tj = (J2 aa * e aa 'b\ ma b) ma ,) and are to 
be determined self-consistently from the mean field type 
Hamiltonian 
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The variational ground state energy, Pmf, of TImf in the 
large-A/" limit is obtained by diagonalising the bosonic 
part of H-mf, integrating over the 2(Af — 1)N S bosonic 
fields &™ ,<T in the action associated with Hmf- One ob- 
tains: 



E 



MF 



AT 



E 

i>j 



n^ijlQijl n^ij \Qij£<7<7' j h-C. 



+ E^(k; Q, A) + E A * (Kl 2 -!-«)• (ii) 

Here w M (k;Q,A) are the positive eigenvalues of the 
bosonic part of HmF) and k = rib/M is kept fixed in the 
limiting procedure 11 ' 12 . The parameter k is a measure 
for the importance of quantum fluctuations: by vary- 
ing k from small to large values, one drives the system 
from the regime dominated by quantum fluctuations to 
the classical regime, i.e., from the disordered into the or- 
dered region. Finally, the GS is obtained by finding the 
saddlcpoint of Emf in the space of the variables Qij and 
x\ subject to the constraints 



dE UF (Q,X)/dXi = 0. 



(12) 



In addition to the GS itself, the spin-spin correlation 
function Gij = (Si ■ Sj) in the GS is an important piece of 
information. In particular, by considering its behaviour 
in the limit \i — j\ — ► oo, one can distinguish between 
LRO and SRO. According to Sachdev 12 , to obtain in 
the Sp(7V) symmetric approach, the SU(2) invariant ex- 
pression Si ■ Sj must be replaced by the Sp(7V) invariant 
expression 



~ J^Wlbfb]^). (13) 



Within the mean field approach, Gy can then be calcu- 
lated straightforwardly. 



B. The anisotropic kagome antiferromagnet 



operators b"' is symmetric under transformations of the 
projective symmetry group (PSG) that is related to the 
symmetry group of the spin Hamiltonian 7Yakaf (Eq. 1) 
(see Rcf. 20). Wc include two translations, a rotation 
by 7r and a mirror axis orthogonal to the preferred direc- 
tion of the exchange constants (J). Thus generalising the 
treatment of Wang and Vishwanath to our model, we find 
eight mean-field states with different symmetries. Seven 
of them have flux in the sense of Ref. 21 in various cells 
of the lattice. Following the arguments in Rcf. 21, we ex- 
clude all flux-carrying states and end up with the solution 
(cf. Fig. 8) Pi. 2 , 3 = Qi.2,3, Qi = Qi, and X b = X a . 
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FIG. 8: (Color online) Arrangement of mean field parameters: 
Qi = Qab, Q2 = Qbc, and Q3 = Q ca denote the intra triangle 
bonds, Pi = Q ba i , P2 = Q c y and P3 = Q ac i denote the inter 
triangle bonds. A a , A;,, and A c are the Lagrange multipliers 
needed to implement the constraints on the sites a, b, and c. 

In order to check the flux-argument in Ref. 21, we have 
cxplicitely studied the solution Pi. 2, 3 = — <9i,2,3 and 
found that it is always of higher energy (For J = 1, this 
agrees with the result of Rcf. 12). 

Thus, the expression Eq. (11) can now be cast into the 
form 



Emf 



1. Choice of mean-field variables 

We wish to apply the procedure described above to the 
AKAF represented by the Hamiltonian (1). To render 
the problem of finding the eigenvalues oj^ in Eq. (10) 
and of optimising Pmf tractable, we have to restrict the 
number of variables Qij and . We do so by demand- 
ing that the mean field Hamiltonian TLmf for the spinon where the condensate is written in diagonalized form and 



J|Qi| 2 + 2|Q 2 | 2 -(2A a + A c )( K + l) 
+^E^( k H 1 + Mk)| 2 ) , 



(14) 



k.// 
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t^j(k) are the three positive solutions of 

detD(w) = 0. (15) 

Here, 

6 <"> = (VA? ui )- m 

with 

A = diag(X a , A c , A a ) , (17) 

/ Q 2 (k) -JQi(-k)\ 
Q= -Qa(-k) Q 3 (k) , (18) 
V JQi(k) -Qg(-k) / 

and Q Q (k) = ig Q (>« k / 2 - e -^ k / 2 ) , a = 1, 2, 3, 

(19) 

£1,2,3, see Fig. 2. 



Technically, we find it convenient to use a polar coordi- 
nate parametrisation for the variables Q\, Q2 and A a , A c : 

Qi = Qcos(a), Q 2 = Qsin(a), (20) 

A = Asin(/3), A c = Acos(/3). (21) 

We perform an optimization with respect to the variables 
Q, A, a, /?, as well as condensate densities x M (k) in accord 
with the following algorithm (J and k are kept fixed). 
i. We fix the angles a, (3 and the amplitude Q, and 
first exploit the stationarity condition for -Emf with re- 
spect to A. It is convenient to write the corresponding 
equation in the following form: 

[2sin(/3) +cos(/3)] (k+1) 

— ■k I d 2 k ^|x M (k)| 2 9 A ^(k) 
" Jb.z. 

= 7i[ d 2 kJ2dA^(k), (22) 
"Jb.z. 



2. Technical details of the numerical extremalisation 

Determination of the ground state of the AKAF in the 
considered approximation has been reduced to minimiza- 
tion of the Eq. (14) with respect to two variables Qi and 
Q2, subject to the Lagrange constraints with respect to 
two parameters A a and A c . Being apparently trivial, the 
optimization procedure turns out to be quite involved 
technically. 

First, we find it crucial to consider at least two differ- 
ent chemical potentials. Other than for the spatially 
isotropic KAF, J = 1, we were not able to find a non- 
trivial solution if we used a single A, A a = At, = A c . If A a 
and A c are different, [A, Q] 7^ 0, the Lagrange multipliers 
enter the expressions for the frequencies tu^ non-trivially, 
other than in the case of a global uniform chemical po- 
tential (cf. Ref. 12). In turn, the Lagrange constraints 
cannot be satisfied semi-analytically, and require a nu- 
merical treatment. Second, we choose to work directly 
in the thermodynamic limit N s — > 00 of the model (14) 
by performing a numerical self-adapting integration over 
the BZ. In this limit, the singularities can be integrated, 
and symmetry breaking is signalled by the appearance of 
a finite value of the condensate amplitude (k) at a cer- 
tain wavevector k = q or <f > which characterises the type of 
magnetic order. We mention here that the extremalisa- 
tion of a mean- field energy of the type of Eq. (14) can also 
be achieved by solving the pertinent stationarity condi- 
tions numerically for finite systems, i.e., for finite N s , 
see e.g., Ref. 22. Then, the type of magnetic order has 
to be detected by calculating the structure factor. Third, 
we see that the Eq. (14) has a minimum with respect to 
the physical bond parameters Q\ and Q2 only after the 
elimination of the chemical potentials. In the full Q — A 
space we face an extremalization problem. 



where D, = 8ir 2 /\/3 is the volume of the unit cell. One 
finds that Q and A enter the Eq. (22) only via the ratio 
£ = A/Q. 

The requirement that the frequencies must be positive, 
w /*(k) > 0, defines a lower limit £min(oi,j3) for £: the 
frequencies w M (k) are positive for £ > £, m in(o:, (3); for 
£ = £,min(ct, (3), the lowest mode lo^ 10 vanishes at some 
point(s) ko in the BZ. When this happens, the corre- 
sponding condensate density x W) (ko) can be put non- 
zero, if this is necessary to satisfy Eq. (22). It is im- 
portant to note that in order to determine the actual 
value of £ TO m(a,/3) (as well as those of Q, a and (3) it 
suffices to only consider Eq. (22) at x M (k) = 0, irrespec- 
tive of whether there is condensate, w AlQ (ko) = 0, or not, 
w M (k) 7^ for all k, fi. 

We solve the Eq. (22) for £ numerically in two steps. 
First, we determine £ m m(a,/?): we decrease £ from large 
positive values until the condition td^ (ko) = signals 
that £ = £mm( a )/3)- Second, we set x M (k) = and at- 
tempt to satisfy Eq. (22) in the interval £ > $ ml „(a„8). 
To this end, we set A = £Q in Eq. (14) and vary £ to de- 
termine the extremum of -Emf (i.e., Eq. (22)). We find 
that the extremum is a maximum. If this maximum oc- 
curs for some £ > £ m i n (<x, (3), then Eq. (22) is satisfied 
with x M (k) = 0. If, however, Emf(o:, f3, £,Q, Q) decreases 
monotonously as we lower £ down to £ = £ m m(a!,/3) , 
then the Eq. (22) cannot be solved with x M (k) = 0. In 
this case, a finite condensate density x Po (ko) 7^ 0, is re- 
quired, in order to "compensate" for too large a value of 
the lhs. of Eq. (22). This fixes both £ = £ m m(a,/3) and 
the value x Mo (k ) (cf. sects III B and IV B of Ref. 12). 
ii. Having determined the value of £, we notice 
that the function Emf{<x, P, A, Q) is quadratic in Q and 
bounded from below, which allows an analytical deter- 
mination of Q as the position of the minimum. 
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Hi. Finally, knowing the values of A and Q, we pro- 
ceed by a numerical extremalization of Emf with respect 
to the angles. The calculations show that Emf as a func- 
tion of the angle (3 possesses a maximum, and a minimum 
as a function of the angle a after f3 has been eliminated. 
Thus, the variational energy Emf is bounded from below 
in the variables Q\ and Q2, as expected. 
iv. We iterate this procedure (i)-(iii) until conver- 
gence is achieved. 



C. Numerical results of the Sp(N) formalism 

The results of the Sp(7V) approach in the large- AT limit 
are summarised in the zero temperature phase diagram 
of the AKAF, Fig. 9. The central part of the phase dia- 
gram is occupied by the incommensurate (IC) phase with 
LRO at sufficiently small 1/k. The phase boundary that 
separates the region with SRO from the region with LRO 
was found by checking whether for a given pair of J and 
1/k the lowest branch of the one spinon spectrum w^(k) 
has zeros in the BZ or not, i.e., whether there will be 
condensate at one or several points in the Brillouin zone 
or not. As one might expect, LRO is maximally sup- 
pressed by quantum fluctuations for J = 1, which is the 
case of maximal frustration. 

For J = 0, the exact quantum ground state of the 
AKAF is ferrimagnetic (FM) according to the Lieb- 
Mattis theorem 15 . In this state, the expectation value 
Qi which measures the singlet weight on the horizontal 
bonds vanishes. As shown in Fig. 10, our Sp(7V) calcu- 
lations recover this exact result and extend it to a finite 
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FIG. 9: (Color online) Phase diagram of Hakaf as obtained 
in the Sp(A/") approach. Symbols and lines, respectively, de- 
note numerical and analytical results for the phase boundaries 
(see text, Subsec. IIIB2 and Appendix B). Quantum fluctu- 
ations increase along the vertical axis. LRO: Long Range 
Order; SRO: short range order; FM: ferrimagnet; IC: incom- 
mensurate phase; DC: decoupled chains. At J = 1, the results 
of Ref. 12 are recovered. Incommensurate order (see Fig. 12) 
occurs between the boundaries of the ferrimagnetic phase ( x ) 
and of the decoupled chain phase (*). 
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FIG. 10: (Color online) Mean field parameter Qi as function 
of the anisotropy. 



interval < J < Jf e m{n), which narrows as 1/k in- 
creases. The parameter Q2, which measures the singlet 
weight on the diagonal bonds, is independent of J in this 
interval; its value decreases as 1/k increases (see Fig. 11). 
Remarkably, the FM state retains its LRO in its entire 
region of existence. 

As J is increased beyond Jf e rri{^), Qi increases in the 
manner of an order parameter at a second order phase 
transition. At the same time, the parameter Q 2 begins 
to decrease, and eventuallly it drops to zero at some 
J = Jdc{ k )- Thus, the large- approach predicts the 
existence of a decoupled-chain phase in the region above 
the phase boundary Jdc{k)- Q2 decreases to zero con- 
tinuously so that the phase transition at Jdc{k) appears 
to be of second order again. 

Both LRO and SRO phases may be characterised by an 
ordering wave vector q or d = 2k m i„, where k m i n is that 
wave vector at which the one-spinon excitation spectrum 
cj M (k) has its minimum. The static spin structure factor 
S(q) develops a peak at q_ or d- In Fig. 12, we display 
the rr-componcnt of the ordering vector q* rd = q* rd {J) 
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FIG. 11: (Color online) Mean field parameter Q2 as function 
of the anisotropy. 
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FIG. 12: (Color online) Ordering wave vector q% ri as function 
of the anisotropy. 



(llrd = 0)- At the kagome point J = 1, |<j£ r(i | = 47r/3 is 
independent of the value of k. For 1/k < 3, the behaviour 
of as a function of J is as expected: as J increases, it 
increases monotonously until the phase boundary Jdc{k) 
is reached and remains constant inside the DC phase. 
However, for 1/k > 3 the function q° rd (J) develops a 
minimum at J f=a 1.5, which becomes more pronounced 
as 1/k increases. 

In Sec. IIIB2 we emphasised that contrary to previ- 
ous applications of the large- M approach to spin models 



11,12,19 



we 



on kagome and anisotropic triangular lattices 
found it essential to consider two chemical potentials A a 
and A c here, one for the spins on the horizontal lattice 
lines (A a ) and one for the middle spins (A c ). We display 
the values of these parameters as functions of J in Fig. 13. 
We have no physical explanation for the behaviour of A a , 
A c as functions of J and k but it is gratifying to see that 
A a = A c at J = 1 independent of k in accordance with 
earlier work . 

As indicated above, along with numerical study of Eq. 11, 
we performed extensive analytical calculations, both to 
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FIG. 13: (Color online) Lagrange multipliers A a , A c (chemical 
potentials) as functions of the anisotropy. 



corroborate the numerics and to obtain additional insight 
into the problem. Details of the analytical techniques are 
presented in Appendix B. Here we state that we were 
able to analytically determine Sp(A/") phase boundaries 
between the SRO and LRO DC phase, between the DC 
and IC phase, and between the FM and IC phase, see 
Fig. 9. Moreover, our analytical calculations allowed us 
to explicitly confirm the existence of LRO inside the FM 
phase and immediately to the right of the FM-IC phase 
boundary. Likewise, the regions with SRO and LRO in- 
side and immediately to the left of IC-DC phase bound- 
ary were determined analytically. This was achieved by 
evaluating in these regions the Sp(A/") generalisation of 
the spin-spin correlation function (S^-S^u) of the model 
defined by expression (13). (u, v = a, 6, c denote the sites 
of the triangular cells i and j of the model, see Fig. 8). 
On the right-hand side of the FM-IC boundary and in- 
side the FM phase, we find for large distances between 



the cells, |r,- 



»1, 



where 



5« 



3 Mk mm )| 2 fir w = a > b 

2N V (X C + X a ) \_ Xn 



(23a) 



(23b) 



w = c 



and u, v = a,b,c denote the sites of the triangu- 
lar cells i and j of the model, see Fig. 8. Here, 
|^3(k, ri m)| 2 /7V v is the condensate density at k m j„ = 
(-tt,0), |x 3 (k mm )| 2 /A^ v = k, sec Eq. (B8). On the 
FM-IC transition line and inside the FM phase, where 
Eqs. (23) are valid, the parameters A a and A c are not 
independent but can be expressed in terms of a single 
parameter 6, see Eqs. (B4), (B12). The sign pattern on 
the right-hand side of Eq. (23) and the ordering wave 
vector q or£ 2 = 2k m i„ = (— 27r,0) are indeed the proper- 
ties one expects to find for the long-distance behaviour 
of the spin-spin correlation function of a ferrimagncti- 
cally ordered state. Since \x3(k m i n )\ 2 /N v remains finite 
for arbitrarily small values of k, the mean- field Sp(A/") 
approach predicts that this order persists in the ex- 
treme quantum limit of our model, 1/k 1. Together 
with Eqs. (23), the fact that the condensate density 
|^3(k m m)| 2 /iVy remains constant inside the FM region, 
see Eq. (B8), implies that the magnetisation of the FM 
phase remains constant up to the FM-IC phase boundary. 
The same behaviour of the magnetisation of a ferrimag- 
netic phase has previously been observed in an exact- 
diagonalisation study of a one-dimensional kagome-like 
antifcrromagnet 23 . At the FM-IC phase boundary, the 
magnetisation becomes spatially modulated with an in- 
commensurate wave vector q or d = 2k m .;„. 
On the left-hand side of the IC-DC boundary and in- 
side the DC phase we find the following large distance 
behaviour of the spin-spin correlation function: 



(Su C ■ Sj. c ) ~ - cos [2 k mi „ (r.; - r,-)] I - 



2q\_ 



-0 



|z 3 (k mm )| 2 + |^3^kmin)lj 

mm J 

for it, i> 7^ c, c . 



-j 2 



(24a) 



(24b) 



Here, qi and A a denote the saddle-point values of these 
variables obtained from Eqs. (B31), (B32). q 2 is a func- 
tion of qi, determined by Eq. (B36) or by Eq. (B44) 
depending on whether 1/k < 1/k s or 1/k > 1/k s 

(2) 

(k s = 0.181, see Fig. 22). w 3 (k TO j„) is the value of the 
second-order expansion coefficient of the lowest spinon 
frequency w 3 (k), cf. Eqs. (B20), (B21c), at its min- 
imum, and 2 k m i„ is the ordering wave vector imme- 
diately to the left on the IC-DC phase boundary and 
inside the DC phase; it is determined by Eq. (B35). 
|x 3 (k m j„)| 2 /iV v = |.T 3 (-k min )| 2 /A v are the condensate 
densities at the wave vectors ±k TO j n . As is shown in Ap- 

(21 

pendix B, (k OT , n ) remains finite for 1/k > l/re s and 



hence |:c 3 (±k m i„)| 2 /A v vanishes. Thus, (Sj, c • Sj 



0, 



i.e., there is no LRO in this region. By contrast, for 
1/k < 1/k s both, |x 3 (±k m j n )| 2 /iV v and u> 3 2) (k mirl ) van- 
ish when the IC-DC phase boundary is approached from 
the left. However their ratio, which determines the spin- 
spin-correlation function, Eqs. (24), remains finite in this 
limit according to Eq. (B42). Thus, for 1/k < 1/k s , 
Eqs. (24) show that while the chain spins S a , Sj, remain 
disordered, there is long-range IC order between the mid- 
dle spins S c along the IC-DC phase boundary and inside 
the DC phase for sufficiently large k. The middle spins 
occupy the sites of a triangular lattice. Remarkably, the 
correlations between these spins predicted by Eqs. (24) 
arc compatible with the spin pattern 
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FIG. 14: (Color online) The kagome lattice as a triangular 
lattice of downward pointing triangles. The coupling strength 
is J on the horizontal bond and unity on the other two bonds. 



IV. BLOCK-SPIN PERTURBATION 
APPROACH 

The basic idea of the block-spin perturbation theory is 
to calculate the states of small clusters of a given lattice 
exactly and to treat the coupling between these clusters 
perturbatively. The basic building blocks of the kagome 
lattice are triangles. Thus it is natural to consider the 
trimeriscd kagome lattice in which the spins on the down- 
ward pointing triangles are assumed to be strongly cou- 
pled whereas the coupling on the bonds of the upward 
pointing triangles are assumed to be weak, see Fig. 14. 
(Clearly, the exchange of the roles of the upward and 
the downward pointing triangles will not affect the fur- 
ther development to be presented in the current section.) 
The Hamiltonian for this trimeriscd model reads 



H(J,j) = H v (J) + 1 H A (J), 0< 7 <1, 



(26) 



S [cos(2k mm r i )e 2; + sin(2k mm r :) )e y ] (25) where H V (J) {Ha(J)) denote those terms in the 

Hamiltonian (Eq. 1) that act on the bonds of the down- 
ward (upward) pointing triangles. We will determine 
approximate GSs of this trimerised model in different 
ranges of J in a perturbation expansion w.r.t. 7. The 
hope is that the results will provide some qualitative 
insight into the GS properties of the non-trimerised 
model H(J, 1) which is our original model Eq. (1). The 
same strategy has previously been applied sucessfully to 
frustrated spin models by several authors 5 ' 24-26 . 



that would obtain if the middle spins S c were classical 
spins coupled by a classical Hciscnberg model with ex- 
change constant J along one lattice direction and cou- 
plings J' along the other two directions with a ratio J' / J 
such that incommensurate order with wave vector 2k m ;„ 
would be established. This persistence of long range or- 
der in the DC phase of the AKAF distinguishes our result 
from the result obtained by Chung et at, Ref. 19, in their 
large-A/" Sp(A/") treatment of the anisotropic triangular 
antifcrromagnet: there the DC phase consists of uncor- 
rected linear spin chains. Qualitative considerations of 
the finite-A/" corrections to the mean-field Sp(A/") result 
led the authors of Ref. 19 to the conclusion that instead 
of the DC phase there is spin-Peierls order in the large- J 
region of their model. In the next section, we will present 
a different approach, a block-spin perturbation theory, to 
get further insight into the properties of the AKAF for 
the physical spin-1/2 case. 



Obviously, the GSs of the unperturbed Hamiltonian 
Tt(J,0) are products of GSs of the individual downward 
pointing triangular plaquettes. The GSs of a single 
plaquette and the corresponding energies are 



1. forJ<l: 
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(|m>-Utt» + (|TIT>-|ITT)) , (27a) 
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FIG. 15: Ground-states of triangular plaquettes. Heavy lines 
depict singlets. The coupling strength is J on the horizontal 
bond and unity on the other two bonds. 



(luT>-itu» + (im>-itu» , (27b) 



£ q ~ a 



-1 + J/4; 



(27c) 



the lattice of these plaquettes which is also triangular, see 
Fig. 14. The zeroth order energy eigenvalues associated 
with the states \A(A1)) do not depend on M: 



Ef {M) = N y {-l + J/A). 



(30) 



Hence, there are in total 2^ degenerate zeroth order 
GSs \A{M)). The single plaquette states |a), \a) satisfy 
the conditions for the validity of the Lieb-Mattis theorem, 
Ref. 15: after a canonical transformation which rotates 
the spins on the sites a and b by tt around the z-axis 
IT) — * i|T)i II) "~ * ~*ll)> an d which leaves the spins on the 
site c fixed the coefficients of all basis states on the right 
sides of Eqs. (27a, 27b) become positive (+1/V6). As a 
consequence, all the GSs \A(M)) satisfy the conditions 
for the validity of the Lieb-Mattis theorem. For J = 
it follows from this theorem that the total magnetisation 
of the exact quantum GS |$ exact) of the Hamiltonian 
Wakaf must be an eigenstate of the total magnetisation 



mtot 



(31) 



2. J > 1: 



w = ^=(\m)-\m)), 



|^> = _(| m )-|UT)) 



e, 3 = £0 = -3/4 J. 



(28a) 

(28b) 
(28c) 



Here, the ket vectors denote the spin state of the 
plaquette in the S z basis. The three arrows inside the 
\cba) symbol denote from left to right the spin direction 
at the sites c, b and a of the plaquettes in Fig. 15. The 
states | a) (|ci:)) and \(3) (|/3)) have total z-spin 1/2 
(—1/2). They can be depicted graphically as shown in 
Fig. 15. From these plaquette states, the zeroth order 
GSs of the Hamiltonian Ti(J 1 7) will be constructed. We 
treat the cases J < 1 and J > 1 separately. 

1. J < 1: Since the states |a), |a) are the GSs of the 
individual downward pointing plaquettes in this case, the 
states 



\A(M))= [] |o,-) n 



a,- 



(29) 



e{M} 



je{N v -M} 



are here the zeroth order GSs of H(J, 7). The set {M} 
is a subset of M out of the A^ v downward pointing tri- 
angles of the 3iV v -site kagome lattice; the subscripts i, j 
denote the position of individual triangular plaquettes in 



with eigenvalue m to t = i.e., \$ exa ct) must be a 

ferrimagnetic state. By continuity, one expects this to 
be the case not only for J = 0, but up to a certain fi- 
nite value of J. This suggests that the state \A(M = 0)), 
c.f. Eq (29), is the appropriate zeroth order GS in this 
case and that the degeneracy of the states | A{M)) is lifted 
by the perturbation 7^a in favour of the state |A(0)). 
To confirm this, we determine the creation energy of a 
flipped plaquette in first order in 7, i.e., the difference 
of the energy of the state with one plaquette spin flipped 
relative to the ferrimagnetic state, and the energy of the 
ferrimagnetic state: 



5^E A (M = 1)=E A{1) 
A simple calculation yields 



E 



HOY 



5^Ea(M = 1) = -7(1 — J) , 



(32) 



(33) 



i.e., to first order, \A(M = 0)), the ferrimagnetic GS is 
stable w.r.t. a flip of a single plaquette spin, as long as 
J < 1. 

As a further check on the stability of the state 
\A(M = 0)), we calculate the dispersion of the excitation 
energy of a propagating single flipped plaquette spin. For 
this purpose, we need to determine the overlap matrix el- 
ements between the state with a flipped plaquette spin 
at the site j and states with a flipped spin at one of the 
neighbouring sites, 



tj,j±s 1 = (a J \(a j± s 1 \~fJS j S j± s 1 \uj)\a :j ±5 1 } = , 



(34a) 



11 



tj,j±S 2 , 3 = ( a j\( a j±62,a\'ySjSj±6 ;i>3 \0ij)\aj±6 a>B ) = --7. 

(34b) 

Here, <5„, v = 1, 2, 3, are the primitive lattice vectors 
of the kagome net, see Fig. 2; they connect the sites of 
the plaquette lattice. Then, by diagonalising the ensuing 
transfer Hamiltonian 

n trans = 7 ]T [\j ( \j + «i)0'l + li - «i)0'l ) 

3 



(li + WI + li-WI) 



~(li + *3><il + li-*3>0'l) }(35) 



where \ j) denotes the state with a flipped plaquette spin 
at site j, we obtain for the kinetic energy of this excita- 
tion: 



*(k) = g7 



J cos(k x ) - cos(%) cos(^^-) 



(36) 



Adding the energy for the creation of a single flipped pla- 
quette spin, we find for the total energy of the excitation 
in the limit of small wave vector k 



w(k) = -7 



(\-j)kl + hl + o(k 4 ) 



(37) 



Obviously, the ferrimagnetic state \A(M = 0)) becomes 
unstable against a propagating flipped plaquette spin 
already at J = 1/4, i.e., much earlier than suggested 
by the excitation energy of a static flipped spin (see 
Eq. (33)). We remark that this bound is independent of 
the actual magnitude of the perturbation parameter 7 
and therefore, the qualitative result may survive in the 
limit 7 = 1. 



with the state with one plaquette spin flipped relative to 
\B(0)). We denote this difference by 5^E B (M = 1) = 
Eb(i) ~ Egmy Surprisingly, we find that the matrix el- 
ements (B(M)\Ha\B(M)} vanish for any choice of M. 

There is no first order correction to the energy E^) M y 

5^Eg(M = 1) = 0. Moreover, we observe that the 
off-diagonal matrix elements (B' (M)\H&\B(M)) , where 
\B'(M)) and \B(M)) contain identical numbers of states 
|/3), |/3) but differ in their distribution over the iV v down- 
ward pointing triangles, also vanish. This implies that, in 
contrast to case (?), a flipped plaquette spin cannot hop 
to a neighbouring site in a first order process. Coupling 
between the spins on the c-sites occurs only in second 
order in 7. It is succinctly described by an effective spin 
Hamiltonian for the c-site spins, which are at the same 
time total spins of the downward pointing plaquettes (see 
Fig. 15): 



H 



= E E { J L S i S i+«v 



i £ c v — 1 



+Ji{SfSf +Su + SfSf + j}.(40) 

Here, Sf , a = x, y, z, denote plaquette spin operators; 
i is the position of a downward pointing plaquette on 
the triangular lattice formed by these plaquettes. The 
exchange couplings jj and Jg~ are given as second order 
matrix elements of 7Ya : 



,11 2 Jv- (B ihi 'i\H A \X) {X\H A \B ihi ^) 

-E 



2 sb - ex 
(B lhl , L \H A \Y) (Y\H A \B lhl , L ) 



2e B 



(41a) 



2. J > I: In this region, the states 

\b(m)) = n ia) n i&> 

ie{M} je{JV v -M} 

with eigenenergy 



(38) 



E 



(o) 

B{M) 



7V v (-3J/4) 



(39) 



are the zeroth order eigenstates of W(J, 7). These states 
consist of free spins on the c-sites and of spin-singlet 
dimers that cover every second bond of the horizontal 
chains of the lattice. We wish to answer the question of 
whether the 2 '-fold degeneracy of these states, which 
results from the degrees of freedom of the free spins, is 
lifted by the perturbation 7H a ; in other words, we want 
to find out whether the middle spins remain decoupled 
from the chain spins. Wc proceed as in case (j). We 
compare in a perturbation expansion w.r.t. 7 the en- 
ergy of the state |-B(0)) with the energy of \B(1)), i.e. 



E 

X 



(B iLi ,i\H±\X) (X\H A \B 



2eb — ex 



(41b) 



and i' 



8 V . Here, the states \Bi, 



are zeroth 



order GSs, Eq. (38), whose spin patterns are identical 
on all sites except for the sites i and i' where the z- 
components of the spins take the values a and a', re- 
spectively; \X) and \Y) are excited states of W v . Of 
course, since the SU(2) symmetry of the original Hamil- 
tonian 7i(J, 7) must be conserved in the derivation of 
7i e ffi the expressions Eqs. (41) must yield identical re- 
sults, Jg = Jg- = Jg v . Non-zero contributions to jjj 
and Jg are obtained if either the same term SjSj' of 
7Ya acts in both matrix elements of the numerators of 
Eq. (41) (two-block contributions) or the terms S^Sfc, 
SfcS,' act in the left and right elements, respectively, 
where the plaquette geometry must be as shown in Fig. 16 
(three-block contributions) . In contrast to the case of the 
isotropic KAF studied by Zhitomirsky 26 , the three-block 



12 



FIG. 16: Configurations of y blocks contributing to the in- 
terblock couplings Jg and Jj^ . Double dashed lines indicate 
that the same term element of Ha acts twice between the V 
blocks at sites i and i' (see also text). 



contributions do not produce three-spin interactions in 
the present case. Rather, they contribute to the ex- 
change interactions Jg and of the Hamiltonian H e //, 
Eq. (40). The evaluation of the expressions (41) yields 



2 1 1 

Sl ~ 1 288 J 



56 



1 



OS 



l 

4. 7 



1 

2J 



(42) 



and 



7_ 
6J 

4J 2 



l 

2,7 



[1 + OiJ- 1 )} 



(43) 



Obviously, these results are useful for J ^> 1. There, 
Tleff represents a spin 1/2 Hciscnbcrg Hamiltonian on 
the triangular lattice of the c-sites with a coupling 
along the 8± direction that is strong (0{^ 2 / J)) in com- 
parison to the couplings in the two other directions 
(C(7 2 /J 2 )). This limiting case of the anisotropic trian- 
gular Hciscnbcrg antiferromagnet (ATHAF) has recently 
been analysed by Starykh and Balcnts with field theo- 
retical methods 14 . These authors find that in the limit 
of strong anisotropy, K = Js 1 /J8 2 ~ * °°! the GS of 
the model Eq. (40) is a collinearly ordered antiferromag- 
net (CAF) as depicted in Fig. 17. Since the ordering 
wave vector Qcaf = (7r, 7t/2) of this phase does not 
evolve continuously from the ordering wave vector q/c 
of the incommensurate (IC) spiral phase of the ATHAF, 




FIG. 18: (Color online) Tentative ground-state of the 
anisotropic kagome antiferromagnet in the limit J ^> 1. Dou- 
ble lines: dimers between the spins on the end points. 



(q/c = {qx(K),0) with -3tt/2 < q x (K) < -tt for 
1/2 < K < oo), they conclude that the IC phase and 
the CAF phase must be separated by a quantum phase 
transition. For the trimcrised anisotropic kagome model, 
Eq. 26, these results have the following implications: 

i ) While in the limit of strong anisotropy J> 1 there is 
long-range collinear antiferromagnetic order among the 
c-site spins, the a- and 6-site spins are paired in singlets, 
see Fig. 18. 

ii) This picture of the GS of the trimerised anisotropic 
kagome model Eq. (26) differs from the result obtained in 
the Sp(A0 approach insofar as for sufficiently large k the 
Sp(A/") approach predicts long range IC order among the 
c-site spins up to arbitrarily large values of J. Thus, if 
the picture of a CAF phase for large anisotropy persists in 
the non-trimerised limit 7i{ J, 7 = 1) of the model Eq. 26, 
one would expect a quantum phase transition between 
the IC phase and the CAF phase of the AKAF similarly 
as for the ATHAF. In closing this section, we remark 
that the calculation that led to the effective Hamiltonian 
tleff, Eq. 40, i.e., to the coupling between the c spins 
in the strongly anisotropic limit, shows clearly that this 
coupling arises from quantum fluctuations of the a and b 
spins. 



SUMMARY AND DISCUSSION 



FIG. 17: (Color online) Collinear antiferromagnetic state 
(CAF) on the triangular lattice 14 . 



In this work, we have studied the ground state (GS) phase 
diagram of the quantum Heisenberg antiferromagnet on 
the kagome lattice with spatially anisotropic exchange 
(AKAF). The model is relevant for a description of mag- 
netic properties of volborthitc, which is a natural reali- 
sation of a spin 1/2 antiferromagnet consisting of weakly 
coupled slightly distorted kagome layers. A small mono- 
clinic distortion along one of the three lattice directions 
causes the exchange coupling along this direction, J, to 



13 



differ from the couplings in the other two directions, J' , 
which we set equal to unity, cf. Fig. 2. We have inves- 
tigated the problem in the full range of the anisotropy, 
< J < oo, using three different approximate methods: 
the classical and semiclassical approach, the mean-field 
Sp(7V) approach, and a block-spin perturbation theory. 
The case J = 1 is the much studied isotropic kagomc 
antifcrromagnet (KAF). Exact diagonalisation studies of 
this model 3,4 are available. Their results speak conclu- 
sively in favour of a spin liquid ground state 1 . This view 
is supported by block-spin approaches 5,6 . Conflicting re- 
sults have been found in Refs. 27-31, where various va- 
lence bond crystal (VBC) states are proposed as ground 
states of the KAF. However, a recent comparison of the 
exact spectrum of the 36-site sample of the KAF against 
the excitation spectra allowed by the symmetries of these 
states, casts doubts on their validity 32 . 
Within the whole anisotropy range, the case J = is spe- 
cial, since it allows for an exact characterisation of the 
quantum GS as ferrimagnetic (FM) with a total mag- 
netisation of M = S N s /3 for a system of N s spins of 
magnitude S. In the classical picture, this state corre- 
sponds to a unique staggered layout of spins with a non- 
zero net magnetisation of the lattice unit cell (cf. Fig. 3). 
In the classical limit, the ferrimagnetic ground state sur- 
vives up to J = 1/2. For J > 1/2, the "chain" spins 
(i.e., spins coupled by J) begin to tilt gradually towards 
the middle (remaining) spins (see Fig. 5). This allows for 
a formation of a large degenerate manifold of canted spin 
states. In contrast to the isotropic case J = 1, where 
the degeneracy grows exponentially with the system size 
N s , its growth is weaker: 2 1 - 26v/ ^ for J ^ 1. This im- 
plies that there must be an increasingly large number of 
classical low energy configurations as J approaches unity. 
In the linear semiclassical approximation, the spin-wave 
spectrum has one zero-frequency mode for each point 
of the magnetic Brillouin zone (BZ). The spectrum is 
identical for the different canted states for all J > 1/2. 
Thus, in this order of the semiclassical approximation, 
no order-by-disorder mechanism appears that would se- 
lect one particular state or a particular group of states 
from the manifold of canted states as true ground states. 
In the limit J — > oo, the frequency spectrum of non-zero 
modes gradually takes the shape of the spectrum that 
one would expect for a set of uncoupled antifcrromag- 
netic spin chains parallel to the strong-J direction. No 
qualitative change from the set of canted spin states to 
the set of decoupled chains at a finite value of J is found. 
We have further explored the nature of the phases at var- 
ious J exploiting the mean field (MF) Sp(jV) approach, 
that incorporates the effect of quantum fluctuations not 
only perturbatively, but self-consistently. The strength 
of quantum fluctuations is controlled by a parameter k, 
which is the analogue of the spin value S in the origi- 
nal SU(2) symmetric model. In fact, for AT = 1, when 
the Sp(l) symmetric model is equivalent to the SU(2) 
model, k = 2S. For general Af, this last identity does 
not hold, but k is still a measure for the importance of 



quantum fluctuations that are strong for k«1 and weak 
for k 3> 1. In the MF Sp(A/") approach, the nature of the 
phases that occur can be read from the values of the 
mean field parameters Q\ and Q2 and from the spec- 
trum of the bosonic spinon excitations. While the mean 
field parameters Q\ and Q2 (cf. Fig. 8) are the GS ex- 
pectation values of singlet bond operators, the structure 
of the spinon spectrum, o> AI (k; Q, A), determines the ex- 
istence or non existence of long-range order (LRO): If 
the spectrum becomes gapless at some wavevector q or d, 
a Bose condensate will form and a modulated structure 
with the wavevector 2q or£ 2 will acquire LRO. 

As was to be expected, the phase diagram of the AKAF 
obtained by the MF Sp(A/") approach contains an incom- 
mensurate (IC) phase in the vicinity of the isotropic point 
,/ = 1 which is ordered for sufficiently large k according 
to this approach, see Fig. 9. Qualitatively, we may gauge 
the value of K against the spin length S by looking at 
the line J = 1 of the phase diagram which is the location 
of the Sp(A/") analogue of the isotropic SU(2) symmet- 
ric kagome model: since, as we have argued above, the 
SU(2) model is disordered for S = 1/2, we may con- 
clude from Fig. 9 that the value of 1/k that corresponds 
to S = 1/2 must be greater than two. Somewhat sur- 
prisingly, the FM phase remains long-range ordered for 
arbitrarily small k. This may reflect the fact that in 
the SU(2) version of the model, the FM phase is ordered 
even for the smallest physical spin value S — 1/2. A new 
feature of the phase diagram is the prediction of a decou- 
pled chain (DC) phase for large enough J, which has no 
classical analogue. In this phase, the chains of strongly 
coupled a- and 6-site spins show no magnetic order. The 
c-site spins which are interspersed between these chains 
and which occupy the sites of triangular sublattice are 
decoupled from the chain spins. However, they may or 
may not exhibit long range order among themselves de- 
pending on the magnitude of k. Remarkably, the spin- 
spin correlations, whose asymptotics were obtained ana- 
lytically, are compatible with the spin-spin correlations 
of an anisotropic classical Heisenberg antiferromagnet on 
the triangular lattice whose exchange couplings differ in 
one direction from those in the other two directions. 

In order to tackle the problem of the GSs of the AKAF 
from a third corner, we have used a block-spin pertur- 
bation theory. This method has the advantage of being 
applicable directly to the spin 1/2 version of the model. 
In applying this approach, one has to initially group the 
spins of the model in clusters. For the kagomc lattice, it 
is natural to choose the spins around either the upward 
or the downward pointing triangles as clusters of strongly 
coupled units and to consider the coupling between these 
clusters, 7, as the small expansion parameter. Thus one 
trimerises the original model (see Fig. 14) and in so do- 
ing, one breaks the translational invariance of the original 
model. In the zeroth order of this expansion, two regions 
can be distinguished by the cigenencrgics of the individ- 
ual trimers: J < 1 and J > 1. For sufficiently small J, 
one recovers the FM state as the GS in first order w.r.t. 7. 
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For J > 1, there are no first order corrections to the en- 
ergy. Following an earlier application of the block-spin 
technique to the isotropic KAF 26 , we determine for J > 1 
in second order in 7 an effective Hamiltonian TL e f f for the 
block-spins which can be identified as the middle spins of 
the original model and that occupy the sites of a trian- 
gular lattice. Ti-eff is a Hcisenbcrg Hamiltonian with a 
coupling J$ 1 of the order of 7 / J along the 8\ direction 
(cf. Fig. 2) and couplings J$ 2 = Js 3 of the order of j 2 / J 2 
along the other two directions. The calculations that lead 
to these results show clearly that the couplings between 
the c-spins of the AKAF are due to fluctuations of the sin- 
glets between the a- and 6-spins into excited states. In a 
very recent field theoretical study, Starykh and Balents 14 
arrive at the conclusion that for Jg 1 ^> Js 2 ,s 3 , the ground 
state of the anisotropic triangular antifcrromagnet rep- 
resented by Ti-eff is a collinearly ordered spin state, see 
Fig. 17. Then, together with the singlet dimers between 
the a- and 6-spins of the downward pointing triangles, the 
state depicted in Fig. 18 emerges as the candidate ground 
state of the AKAF in the limit of large anisotropy J> 1: 
while nearest neighbour spins on the strongly coupled a-b 
chains form singlets and decouple magnetically from the 
spins on the c sites, the latter order in a collincar antifcr- 
romagnetic structure. This structure cannot be obtained 
by a continuous deformation of the spiral IC structure 
that is predicted by the Sp(7V) approach and is believed 
to prevail for sufficiently large k in the region of moder- 
ate anisotropy. As a consequence of the trimerisation, the 
state depicted in Fig. 18 breaks the translational symme- 
try of our original model, Eq. (1). If this state survives as 
the ground state of the non-trimerised model, i.e., when 
the expansion parameter 7 approaches unity, then, ow- 
ing to their incompatible symmetries, the spiral IC phase 
and the large J phase of our model must be separated by 
a quantum phase transition. 
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APPENDIX A: GROUND STATE DEGENERACY 
FOR GENERAL J 

We first derive the constraint on the chiralitics that leads 
to the reduction in the number of degenerate ground 
states for general J relative to the special case J = 1. 
Let xij ' • • 1 X6 be the chiralities of the six triangles sur- 
rounding one of the hexagons of the kagome lattice, and 
let 0i , ■ ■ • , 6 denote the angles that define the direc- 
tions of the spin vectors on the six corners of the hexagon, 



see Fig. 19. 

Then, as is seen in Fig. 19, the following relations between 
the angle 0i , and the angles 02 • • • , 06 are an immediate 
consequence of these definitions: 



'2=01 - 


9X2, 






(Ala) 


'3 = 01 - 




-X3), 




(Alb) 


'4 = 01 - 


Hx2- 


~X3)^(2n 




(Ale) 


'5 = 01 ~ 


Hx2- 


- X3 + XS) " 


-(27T-20) X 4, 


(Aid) 


'6 = 01 - 


Hx2- 


" X3 + X5 4 


■X6)-(2tt-20) X4, 


(Ale) 


and 


06 = 


0i + (2tt - 


20)X1 ■ 


(Alf) 



From the last two of these relations it follows that the 
chiralitics xi, ■ • • , X6 arc constrained by the sum rule 

X2 + X3 + X5+X6- 2xi -2X4 = 0. (A2) 

For the isotropic kagome system, J = 1, 6 = 2ir/3, one 
finds instead of the constraint (A2) the sum rule 

6 

y~^Xj = n where n = 0, 1, 2 (A3) 

which is obviously less restrictive than (A2). 

Next, we present the arguments that lead to the estimate 

7Vg™s°(iV v ) < 2 a V^ with a < 3 (A4) 

for the number N^f so (N^) of classical GSs of an 
anisotropic kagome AF with N-^ downward pointing tri- 
angles (the number of sites is 3A^ V ). Any planar config- 
uration of a cell of the kagome lattice can be constructed 
by decorating the successive rows of up and down point- 
ing triangles with chirality values x = il starting with 
the first row. We consider only square cells with y/N^ 
rows with y^N v downward pointing triangles. Then, 
each row consists of 2y / 7V v triangles, see Fig. 20. 

Obviously, there are 2 2 V Jv v wa y S to decorate the first 
row. Disregarding certain exceptions, which will be dis- 
cussed below, one can, for a given configuration of the 
first row, choose the chirality of an arbitrary triangle 
of the second row to be either +1 or — 1. After this 
choice has been made, the constraint (A2) fixes the chi- 
ralitics of all the remaining triangles of the second row 




FIG. 19: Chiralities around heaxagonal plaquette 
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V V 



FIG. 20: Example of a chirality distribution; dark and light 
shaded triangles represent positive and negative chirality, re- 
spectively. Chirality configurations in boxes fix the chirality 
distribution of the row above them uniquely. An empty cir- 
cle inside a triangle indicates that its chirality can be chosen 
freely to be positive or negative. 



uniquely. Proceeding in this manner from row to row 

one would generate 2 2 "v/ A ' v ■ 2"^ Nv distributions of chi- 
ralitics over the downward pointing triangles of the 
cell. For finite lattice cells, the requirement of periodic 
boundary conditions imposes further constraints on the 
number of possible chirality distributions in these cells, 
but the effect of these constraints will become negligible 
in the thermodynamic limit iV v — > oo. However, there 
is a further reduction of the number of possible chirality 
distributions: For a given distribution in a row it is not 
always possible to find two distributions for the successive 
row which both satisfy the constraint (A2). If in a row 
the lower half of a hexagon of the next row is decorated 

by chiralities in the manner — H or H h (see boxes 

in Fig. 20), then the chiralities of the next row are fixed 
uniquely. This reduces the number of possible chirality 
distributions. Obviously, this reduction of the number 
of possible chirality distributions survives in the thermo- 
dynamic limit so that the exponent in (A4) is less than 
3 ^/A r v , the value one would have expected without this 
reduction. We have calculated the number of distribu- 
tions for cells of up to N v = 13 x 13 and have found the 
value a ~ 2.18 for the constant in the expression (A4), 
see Fig. 21. 

As we have mentioned above, the sum rule (A3) which 
applies for the isotropic kagome AF is less restrictive than 
the sum rule (A2). Consequently, the number of chirality 
distributions in the isotropic model 33 , 



N}ls ~ 1.1833 3JV "= 



(A5) 



is larger than in the anisotropic model. Since the tran- 
sition from the anisotropic model to the isotropic model 
happens through a continuous variation of the coupling 
constant J, there should be a continuous transition be- 
tween the numbers of GS configurations in these two 
cases. Presumably, this transition implies that the den- 
sity of low-energy states of the anisotropic model in- 
creases exponentially with an exponent ~ */ iV v so that 
for J — > 1 a sufficient number of states collapses to the 
GS to bring about the transition between the laws (A4) 
and (A5). 




FIG. 21: (Color online) Number of chirality distributions, 
Nas s °, of cells of up to iV v = 13 x 13. Dotted lme:_jnin = 
2 W iV v In 2 (lower bound); dashed line: max = 
(upper bound); full line: ln(7Vg5 ISO ) = 0.65 + 2.18 
(linear fit to the numerical results). 



APPENDIX B: PHASE BOUNDARIES 

The FM phase and the DC phase are chacterised by 
the vanishing of the parameters Q\ and Q2, respec- 
tively. Our numerical results in section III C show that 
at the respective phase boundaries, Q\ and Q2 decrease 
to zero like order parameters at second order phase tran- 
sitions. This suggests that we expand the mean field 
energy Euf, Eq. (14), w.r.t. either Q\ or Q2 in the 
manner of a Landau- Ginzburg (LG) expansion and de- 
termine the phase boundaries and the properties of the 
FM and the DC phase from this expansion. We write 
(«), 



E M f/{N v N) = egi(Qa) where 

4 Q i(Q Q ) = e a +r a \Q a \ 2 + g a \Q a \ 4 + 0(\(Q a \ 6 ) . (Bl) 

The coefficients e a , r a and g a are functions of the vari- 
ables k and J, of the parameters A a , A c and of Q2, Qi 

for a = 1, 2, respectively. The saddle point of e£g(Q a ) 
w. r. t. A a , A c and Qp, [3 ^ a, determines the physical 

values of these parameters. For e£g(Q a ) to qualify as 
a bona fide Landau- Ginzburg energy describing a second 
order phase transition with Q a playing the role of an or- 
der parameter, the coefficients g a have to be positive at 
the saddle point. For g\, i.e., inside and on the bound- 
ary of the FM phase, this follows from the numerical 
result: Q\ is found to remain zero for all J < Jf(k). By 
contrast, we have no numerical results for J > Jdc( k )i 
i.e., inside and on the boundary of the DC phase. There- 
fore, we need to show by analytic means that 172 > 0. 



1. The FM phase and the FM-IC phase boundary 

Since, as we have just remarked, we know that g\ > 0, 
the remaining task is to determine the coefficients e\ and 
ri of e^Q. To this end, we have to expand the mean 
field energy Emf, Eq. (14), w.r.t. Q\ which amounts to 
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expanding the frequencies w M (k) w.r.t. Q\. As can be in- 
ferred from the expressions (16), (18) the frequencies de- 
pend on Qi only through the combination e 2 = J 2 \Qi\ 2 - 
Therefore, we write the expansion in the form 

^(k;e)=c(°)(k)+e 2 4 1 )(k) + 0(e 4 ) 



(0) 



(k) 



(k)+0( £ - 4 ) (B2) 



,-,(o) 



■ J n / A + 



(k) 



with A + = (A a + A c )/2, 
\+d e 2 w A1 (k;e)| e=0 and e = e/A + . 

Here, the introduction of the "dimcnsionlcss" quantities 
_ (i) _ 

ujji , e looks like an unncccessary complication but it will 
help to keep expressions further below simple. Setting 
Qi = in the matrix D(w), Eq. (16), and solving Eq. (15) 
for uj we find 



Ld[ 0) (k) =W F (k) + 5, 
^ 0) (k) = 1-6, 
4 0) (k) = w F (k)-S. 



Here 



5 = X-/X+ with A_ = (A c - A Q )/2 



(B3a) 
(B3b) 
(B3c) 



(B4) 



and 



u*(k) = yT-^ 2 [sin 2 ( S2 /2)+sin 2 ( S3 /2)] (B5) 



with 



?2 = IQ2I/A+ 

and s a = S a k, 



2, 3 (sec Fig. 8) . 



(B6) 



From our numerical results, Fig. 13, we know that A c > 
A a and hence 6 > 0. Therefore, 0)3 (k) < u)i° 2 (k), 
and hence, if the minimum of lj^ (k) vanishes at the 
point kmin in the Brillouin zone, Wi 2 (k TO j n ) will be fi- 
nite. Thus, since condensate can only occur when one of 

the frequencies u>j? , fx = 1, 2, 3 vanishes there may be 
a finite condensate density |x3(k m i„)| 2 , but the densities 
\x±\ and |x 2 | will certainly be zero. With these remarks 
and with the above results for uijj 3 we find from Eq. (14) 



following three conditions to be fulfilled: 
d ei /d\- = : 



dei/d\+ = 



k: 



(B8) 



2A+g 2 2 - -K - 1 + E 2 (g 2 ) + ^w F (k min ) = 0; 



(B9) 



dejdq 2 = : 

2 A+ - [K 2 (q 2 ) - E 2 (q 2 )} 

Q2 w F (k min ) 



= 0; 
(BIO) 



with K 2 (<f 2 ) = - / ds 2 - ds 3 Wf(k)" 1 , 
T Jo T Jo 

1 r 1 r 

E 2 (g 2 ) = - / ds 2 - ds 3 w F (k) . 
n Jo 7T Jo 

(Bll) 

According to Eq. (B8), condensate must be present in 

the FM region. This requires that lu^ (k m in) vanishes. 
From Eq. (B3c) it is seen that k min = (— tt, 0), so that 

,-,(0) 



(kmin) — 0, if 



w F (k min ) = \J\-2q 2 



A_ 
Al 



(B12) 



Within the FM region and on the FM-IC boundary 
(i.e., for Qi = 0) the saddle-point values of q 2 , A+, 
and A_ are then determined as functions of k by the 
Eqs. (B9), (B10) and (B12). Remarkably, within this re- 
gion these quantities are independent of the value of the 
exchange constant J. The solution of these equations 
shows that < q~ 2 < 2/3 for < k < 00, cf. Figs. 11, 13. 
The FM-IC phase boundary is the solution of n (k, J) = 
(cf. Eq. (Bl), where 



n = de^l/dQ 2 



(B13) 



with e^l (E M f) horn Eq. (14). 
We obtain 



ei/A + = 2A+<7 2 -(3-(5)(K + l) 



+ 



rain ) |^3 \K-rain 



,(o) 



(k) 



-,(0) 



(k) 



(B7) 



Stationarity of e\ w.r.t. A_, A+, and q 2 (which is equiv- 
alent to stationarity w.r.t. A a , A c ), and Q\ requires the 



J 2 



1 



1 1 



+ - — lim 

A+ k^k mi 



S2 + S3 



ftW(k) 



5m 2(f2+f3 , ^V (k) ) (B14) 



with n« (k) = -w^ (k) - 4 1} (k) 



-,(!) 



(k). 



To obtain the expansion coefficients lu^ (k) which appear 
in the last equation, we solve Eq. (15) to first order in 
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the expansion w.r.t. e 2 . We find 



a. Expansion of e^l (Q2 



nM(k)= 



w F (k)(w F {k) + 1 - 28) 



q 2 A sin 2 (sa.) sin 2 ls 



w F (k) + 1 
, (2w F (k) + 1-8) 



(w F {k) + 1)(1 - 8) (w F {k) 2 -S 2 ) 



(B15) 



and 



lim 

k^k mi 



sin 



S2 + S3 



-,(!)/ 



1 1 



28 1 



(B16) 

With these results Eq. (B14) can, in the thermodynamic 
limit, be cast into the form 



ri_ 

J 2 ~ 
where 

^3(92) 



1 

J 



^3(32 



K 

Al 28 1 



1 1 + 8 



1 
J 



1 



J f (k 



dsn 



ds-i 2 sin 



(B18) 

Then, with (72 = 92 (k) and A± = A±(/t) as obtained from 
Eqs. (B9), (BIO) and (B12), the condition n = is 
an equation for the FM-IC phase boundary J = J f (k) 
which yields the graph shown in Fig. 9. As we have 
mentioned above, inside the FM region, i.e., for J < 
J f (k), the saddle-point values of the quantities q and A± 
and hence of Q 2 , A a , A c and |a?3 (kmin ) | arc independent 
of the exchange coupling J, i.e., they retain the values 
they attain on the FM-IC phase boundary, cf. Figs. 11, 
13. 



2. The DC phase and the IC-DC phase boundary 

Proceeding in exact analogy to the development in the 
previous subsection we now expand Emf I (-^y N) in 
powers of |Q2| 2 - However, instead of working with the 
variables Q±, Q 2 , A a , A c we work with qx, Q 2 , A a , q 2 here, 
where 



92= 



J\Qi\ 

Xa 
IQ3] 



VA a A c 



(B19a) 
(B19b) 



The replacement of \Qi\ is purely a matter of conve- 
nience. By contrast, the replacement of variables Q 2 , A c , 
which according to the numerics vanish simultaneously 
as J approaches the IC-DC phase boundary, by the pair 
Qii qi leaves us with only one vanishing variable, since, 
as will be seen below, q 2 remains finite throughout. 



We write 



u, M (k) = W W(k) 



^ 2) (k)g 2 



(k)Qj +£){Ql 



(B20) 

and determine the coefficients to^ , n = 1, . . . 4, by solv- 
ing Eq. (15) for uj iteratively. We obtain 



wi 0) (k)+w^(k) = 2X a w DC (k). 



,(°) 



,(°) 



<[ 2) (k)+ W ^(k) = -f- 



( 2 )/^\ - 



1 1 — cos k x cos k y 



X a w DC {k) 



4 2 \k) - 



[C(k) 2 - D(k) 2 } 



1/2 



, (B17) Here ' 



C(k) = l-q 2 2 



1 — cosk x cosk v 
w DC (k) 2 



D(k) = qi smk x q 2 



cos k x — cos k v 
w DC (k) 2 



0, (B21a) 
(B21b) 

(B21c) 

(B22) 
(B23) 



w DC (k) = \Jl - q\ sin 2 k x 



(B24) 



The coefficients u)jx'(k), \i = 1, 2 ,3, will only be needed 

in the determination of the coefficient g 2 of the fourth 

( 2) 

order term of e K L ^(Q 2 ) which will be discussed later. We 

will first concentrate on the determination of the zeroth 

order term, e 2 , and of the coefficient r 2 of the second 
(2) 

order term of c L q(Q 2 ). Under the assumption that g 2 is 
positive, this will provide us with an expression for the 
IC-DC phase boundary. 

With the above expressions for and wl , 

u = 0, 1, we obtain for the coefficients of the Landau 
Ginzburg energy from Eqs. (14), (Bl) 



e 2 {qi,X a ) = 



\ 2 n 2 
A a ffl 

J 



2X n 



1 + K- — VfflDc(k) , 

ivv k 

(B25) 



r 2 (qi, 92, A a , |x 3 (k mm )| 2 ) 

11/ 

-V q 2 



£ 1 

k 



A a ^ w DC (k) 



1 



Y,^ ] (k) + \x 3 (k mm )\ 2 J 3 2 \k 



(B26) 
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These are valid for arbitrary values of the parameters q\ , 
A a , q 2 j and |a;3 (k miT j) | . In the next subsection, we will 
calculate their saddle point values for given Qi and thus 
fix the parameters. Here, we have only allowed for the 
existence of a condensate component |x3(k m .;„)| 2 . This 
is justified since, as Eqs. (B20) and (B2) show, uj^ < wi,2 
for sufficiently small Q 2 so that conceivably W3(k) may 
vanish at some point k m i n in the the Brillouin zone, while 
o>i(k) and oj 2 (k) remain finite at k,„i n , and hence a finite 
condensate density |a; 3 (k mi „)| 2 may occur at this point. 



which determine the saddle point values qf 1 and Aa '. 
( K and E are the elliptic integrals of the first and the 
second kind.) 

(2) 

Next we seek the extremum of e LG w.r.t. q^. Since ei is 
independent of qi we, neglecting terms of order have 



0=d q2 r® 



4. X a 



re + 1 — h(qi, 92) 



b. Saddle point, phase boundary 



1 C(k min ) \x 3 (k 



(B33) 



( 2") 

Next we need to determine the saddle point of e LG (Q2) 
in the space of the variables qx, A a , 52, and |x3(k m i n )|. 
First, the saddle point values of q\ and A a are obtained 
as expansions in powers of Q2, 



X a =X^ + X^Ql+0(Ql), 
0(Q$), 



(0) , (2) n 2 , 
1l V2 



<h=<h 



where , q^ are the solutions of 



(B27a) 
(B27b) 



and 
h(qi, 92) 



2 W 2 1 

- / dk x - I dky ■ 



C(k) 



[c(k) 2 - D(ky 



(B34) 

(In these expressions and in the sequel, we use an abbre- 
viated notation: A a , qi and q2 denote the zeroth order 
quantities A« , gf^ and q^ ■) k TO , n is the location of the 
minimum of (k), 



9A o e 2 =0, 
d qi e 2 =0. 



(B28a) 
(B28b) 



Since the first derivatives of e2 vanish at A a = X^ , q\ = 
qf\ Eqs. (B28), we have 



k v - n ■ 



v 2 ; 



9i 2 



(B35) 

From (B21c) and (B35) it follows that w 3 (k m in) = 0, 
if 



e 2 = 4 J + 4^2 +o(Ql), 



and 



(B29) 



ql = (1 - g 2 )/2. 



(B36) 



As a function of 52 the integral Ii(q\, 52) increases 
monotonously, 



'2 



r { 2 0) + r« QI. 



■O(Qt) 



(B30) 



Here. 



4 and ^2°'' 



are the expressions (B25) and (B26) 



with A a and qi replaced by X^ and q\"' . The fourth 

order term of e2, Eq. (B29), and the second order term 

(2) 

of T2 contribute only to the fourth order term of e LG 
which will be determined later. Therefore, we postpone 
the presentation of explicit expressions for A a , q\^ and 
the ensuing expressions for e 2 and r 2 until later. With 
e2 from Eq. (B25), Eqs. (B28) yield the equations 



re = -K(q^)-l, 

IT 



\(0) 



( 3i 0) ) 



k(^) - E( ff w; 



(B31) 



(B32) 



1 = 0) <h(qi, 52) </i(gi) 



for < q 2 < yj(l - g 2 )/2 
We have defined 



(B37) 



max/i(gi, g 2 ) = h <7?)/ 2 

(B38) 

As is seen in Fig. 22, the graphs of the functions 
re = re(gi), Eq. (B31), and of ii = h(qi) intersect at 
q ls ~ 0.708, re s ~ 0.181. Therefore, in solving Eq. (B33) 
for g2, two cases have to be considered separately: 

i. qi > qi s , re > re s . In this solution ex- 

ists only, if the last term in parentheses in Eq. (B33) is 

( 2") 

positive. This requires that oj^ (k 
has been discussed before, |x3(k TO ,, 

tio |x3(k m i„)| 2 /a;3 2 (k m i„) would vanish otherwise. The 



m ) =0 because, as 
and hence the ra- 
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FIG. 22: (Color online) Ii and k + 1 as functions of qi. 



condition Lu 3 (k mi n) = implies that gj = (1 — qf)/2, cf. 
Eq. (B36). Using this result and Eq. (B32) to eliminate 
q 2 and A a from Eq. (B26) we find 



(see Eq. (B38)). Consequently gf < (1 - q\)/2 so that 

(2) 



[X3 



3 (k m i„) > and hence no condensate can develop, 
2 = 0. Then, Eq. (B33) yields the equation 



h(qi, 92) 



1 



(B44) 



which replaces Eq. (B36) and determines q 2 as a function 
of qi, q2 = Q2 ) ■ Then, proceeding as in case (?) one 
finds for the IC-DC phase boundary in the region k < n s 



JDc(K)=[(K+l)(l + qa) + h(qi, 92) 
ql 4 [K( 9l ) - E( 9l )] • 



(B45) 



.(0) 



2 1 



Jdc(k) 
J 



(B39) 



where 



Jdc(k) 



with 



( K + l)(3-q()/2 + I 2 ( qi , x /(l-qj)/2) 
2 



q 2 n 



(l-ql) 4[K( 9l )-E( (?1 )] 



(B40) 



^2(91,92) 



-7T 2 , 

^ -10 



dk x / dk y [C(k) 2 - D(k) 2 ] 
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(B41) 

is the IC-DC phase boundary for k > k s , i.e., in the 

region where the ratio \xz(\t m in)\ 2 1 (\L m in) is finite. 
According to the discussion at the end of section III C, 
cf Eq. (24) , this is the region where LRO prevails along 
the decoupled chains, cf. Fig. 9. 

In the development leading to Eq. (B40) for the phase 
boundary, we have not needed the solution of Eq. (B40) 
explicitly, but we note it here for completeness: 



1 C(k 

mi n ) |z3(k 



2 i 

^31 



^q\ 4 2 \k mm ) V " 2\ 



1 _ 2 
«f 1+39? 



--l + K-h(qi) > 0. (B42) 



These relations show that while |a;3(k ITt ,- rl )| = 0, the ratio 
|£3(k m i„)| 2 /c4 2) (k min ) remains finite. 

ii. q\ < qi s , K < k s . In this case, we must have 



(B43) 



Here, q\ = qi(n) from Eq. (B31) and q 2 = q 2 (n) from 
Eq. (B44)(with 9l = ft (ft)). 

We note here that inside the DC phase, i.e., for J > 
JdC'(k), where Q 2 = X c = 0, the saddle-point values 
of ft and \ a / J and hence of Q\ are independent of J, 
cf. Eqs. (B31), (B32). Hence the graphs of Qi and A a for 
J < J dc and for J > Juc join smoothly at J = Jdc, cf. 
Figs. 10, 13. Furthermore, it follows from Eq. (B42), that 

the ratio (|x 3 (±k mi „)| 2 /7V v )/(A a gfw 3 2) (k mm )), which 
occurs in the amplitude of the spin-spin correlation func- 
tion, cf. Eq. (24a), is also independent of J inside the DC 
phase and retains the value that it has attained at the 
IC-DC transition line. 



c. Stability of the phase boundary 



In deriving the phase boundary from the condition r\ 
we have tacitly assumed that the coefficient g 2 of the 
fourth order term in the LG expansion, Eq. (Bl), is pos- 
itive. In the remaining part of this appendix we will 
sketch the steps which lead to the conclusion that this is 
indeed the case. 

( 2 ) 

Expanding in the expression (Bl) for e LG the coefficients 
e 2 and r 2 w.r.t. the second order contributions to ft and 

(2) (2) 

X a , ft and Aa , cf. Eqs. (B27) we obtain 



(2) 



where 



rf Q 



+ (92 + g' 2 ) Q\ + 0{Qt) , (B46) 



.92 



EG 



,(4) 



< 4 > 



(B47) 



(2) 

is the contribution to the fourth order term of e L Q that 

arises from the fourth order terms of the frequencies w M 

f 2) 

in the sum in Eq. (14) whereas the contribution to e LG 
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of the expansion of e 2 and r 2 



is 



92=2 



:) x{2) d^e 2 \ d qi d Xa e 2 \ 




d Xa d qi e 2 \ dl a e 2 \ 

d m r 2\ \ 
9\ a r 2 \ ) 




(B48) 



(In Eq. (B48) the notations c^^Iq etc. indicate that af- 
ter the derivatives have been taken the variables qi, X a 

etc. have to be replaced by their zeroth order values 9^ , 

Ai 0) etc.) 

The evaluation of the contribution (B47) is straight- 

foward: the coefficients lo\^\ /i = 1, 2, 3, were obtained 
by solving Eq. (15) for u> iteratively to fourth order. As 
the explicit expressions are rather lengthy and contain 
no direct information, we refrain from presenting them 
here. The sum over k that is required in Eq. (B47) was 
done numerically. g 2 was obtained in the form 



Inserting these results into Eq. (B48) one finds 



9' 2 = -\{d qi r 2 \ ^r 2 \ )M^A 



(B53) 



While the second derivatives of e 2 are obtained straight- 
forwardly from Eq. (B25) the derivatives d qi r 2 \ and 
d\ a r 2 \ have to be calculated separately for the region 
qi < qi s , where there is no condensate, |x3 (kmin.) | 2 = 0, 
and for the region q\ > qi s , where |x3(k m i„)| 2 > 0. Fi- 
nally, the result for g' can be cast into the form 



where 



92 = (*? x\) M' ( 



(B54) 



.92 = -^92{qi) 



(B49) 



where g 2 {q\) is a function of q\ alone which is always 
positive so that g 2 > throughout. Remarkably, no 
explicit dependence on the coupling constant J appears 
in these results. 

The evaluation of g' , Eq. (B48) requires the knowledge of 

(2) (2) 

explicit expressions for qf and Aq . These are obtained 

(2) ( 2 ) 

by expanding e 2 to first order in q\ and X a , inserting 

(2) 

the results into the expression (Bl) for e LG and requiring 
that the terms of order Q\ satisfy the extremum condi- 
tions w.r.t. q\ and A a : 



M' 



with 



1 



A 



-A 



Z?A(«+1-A) \ _ A {2 _ q 2 )A _ K _ 1 



(B55) 



1 2 

A = —x — [Kfai) - E(gi)l (K, E : elliptic integrals) 
qf 7T 

(B56) 

and 



= q[ 2) d 2 qi e 2 \ + X^d qi d K e 2 \ B + d qi r 2 \ ,(B50a) 
= q[ 2) d qi d Xa e 2 l+X^d 2 Xa e 2 \ +d Xa r 2 \ (B50b) 
The solution of these equations reads 




with 



M~ 



r ,5 9l r 2 | 
,<?A a r 2 | 



d qi e 2 \o d qi d Xa e 2 \ 
,d Xa d qi e 2 \ Q d\ e 2 | 



(B51) 



(B52) 



1 - q z 

X q = gf(A-K-l)H 2^-qid qi I 2 (qi,q 2 ) 

?2=g2(9l) 

-e(<n- qis)^^ (K+l-hfa)) , 
1 - qf 

x\ = {\ + 1 )( K + !) - \h{qi,q2, 
qi 92 

92=92(91 ) 



(B57) 



Here O is the step function; the integrals Ii(qi) and 
-^2(91,92) have been defined above, cf. Eqs. (B38) and 
(B41), respectively. After numerical evaluation of these 
integrals, we find that g' 2 = 9 2 {qi) is positive for all values 
of qi. 



1 G. Misguich and C. Lhuillier, Frustration in Two- 
Dimensional Quantum Antiferromagnets (World Scien- 
tific, Singapore, 2004), chap. 5. 

2 R. Moessner, Can. J. Phys. 79, 1283 (2001). 

3 P. Lecheminant, B. Bernu, C. Lhuillier, L. Pierre, and 



P. Sindzingre, Phys. Rev. B 56, 2521 (1997). 

4 C. Waldtmann, H.-U. Everts, B. Bernu, C. Lhuillier, 
P. Sindzingre, P. Lecheminant, and L. Pierre, Eur. Phys. 
J. B 2, 501 (1998). 

5 F. Mila, Phys. Rev. Lett. 81, 2356 (1998). 



21 



6 M. Mambrini and F. Mila, Eur. Phys. J. B 17, 651 (2000). 

7 Z. Hiroi, M. Hanawa, N. Kobayashi, M. Nohara, H. Takagi, 
Y. Kato, and M. Takigawa, J. Phys. Soc. Jpn. 70, 3377 
(2001). 

8 A. Fukaya, Y. Fudamoto, I. M. Gat, T. Ito, M. I. Larkin, 
A. T. Savici, Y. J. Uemura, P. P. Kyriaku, G. M. Luke, 
M. T. Rovers, et al, Phys. Rev. Lett. 91, 207603 (2003). 

9 F. Bert, D. Bono, P. Mendels, J.-C. Trombe, P. Millet, 
A. Amato, C. Baines, and A. Hillier, J. Phys. Condens. 
Matter 16, 9829 (2004). 

10 F. Bert, D. Bono, P. Mendels, F. Ladieu, F. Due, J.-C. 
Trombe, and P. Millet, Phys. Rev. Lett. 95, 087203 (2005). 

11 S. Sachdev and N. Read, Int. J. Mod. Phys. B 5, 219 
(1991). 

12 S. Sachdev, Phys. Rev. B 45, 12 377 (1992). 

13 W. Apel, T. Yavors'kii, and H.-U. Everts, J. Phys. Con- 
dens. Matter 19, 145255 (2007); W. Apel, T. Yavors'kii, 
and H.-U. Everts, ibid., Corrigendum, in press. 

14 O. A. Starykh and L. Balents, Phys. Rev. Lett. 98, 077205 
(2007). 

15 E. Lieb and D. Mattis, J. Math. Phys. 3, 749 (1962). 

16 S. Brehmer, H.-J. Mikeska, and S. Yamamoto, J. Phys. 
Condens. Matter 9, 3921 (1997). 

17 A. B. Harris, C. Kallin, and A. J. Berlinsky, Phys. Rev. B 
45, 2899 (1992). 

18 C. H. Chung, J. B. Marston, and S. Sachdev, Phys. Rev. 
B 64, 134407 (2001). 



19 C. H. Chung, J. B. Marston, and R. H. McKenzie, J. Phys. 
Condens. Matter 13, 5159 (2001). 

20 F. Wang and A. Vishwanath, Phys. Rev. B 74, 174423 
(2006). 

21 O. Tchernyshyov, R. Moessner, and S. L. Sondhi, Euro- 
phys. Lett. 73, 278 (2006). 

22 G. Misguich, B. Bernu, and C. Lhuillier, J. Low Temp. 
Phys. 110, 327 (1998). 

23 C. Waldtmann, H. Kreutzmann, U. Schollwock, 
K. Maisinger, and H.-U. Everts, Phys. Rev. B 62, 
9472 (2000). 

24 V. Subrahmanyam, Phys. Rev. B 52, 1133 (1995). 

25 C. Raghu, I. Rudra, S. Ramasesha, and D. Sen, Phys. Rev. 
B 62, 9484 (2000). 

26 M. E. Zhitomirsky, Phys. Rev. B 71, 214413 (2005). 

27 P. Nikolic and T. Senthil, Phys. Rev. B 68, 214415 (2003). 

28 J. B. Marston and C. Zeng, J. Appl. Phys. 69, 5962 (1991). 

29 C. Zeng and V. Elser, Phys. Rev. B 51, 8318 (1995). 

30 A. V. Syromyatnikov and S. V. Maleyev, Phys. Rev. B 66, 
132408 (2002). 

31 R. Budnik and A. Auerbach, Phys. Rev. Lett. 93, 187205 
(2004). 

32 G. Misguich and P. Sindzingre, J. Phys. Condens. Matter 
19, 145202 (2007). 

33 R. J. Baxter, J. Math. Phys. 11, 784 (1970). 



